Quantum state engineering by a coherent superposition of photon subtraction and 

addition 
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We study a coherent superposition ta + ra) of field annihilation and creation operator acting on 
continuous variable systems and propose its application for quantum state engineering. Specifically, 
it is investigated how the superposed operation transforms a classical state to a nonclassical one, 
together with emerging nonclassical effects. We also propose an experimental scheme to implement 
this elementary coherent operation and discuss its usefulness to produce an arbitrary superposition 
of number states involving up to two photons. 
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I. INTRODUCTION 

Manipulation of light field at the single-photon level 
provides a crucial basis for many important applications 
in quantum information science. In particular, two ele- 
mentary operations on a single-mode field, i.e. photon 
subtraction and addition represented by bosonic annihi- 
lation and creation operators a and a\ respectively, can 
be employed to transform a field state to a desired one [l| . 
For example, the photon subtraction transforms a Gaus- 
sian entangled state (two-mode squeezed state) to a non- 
Gaussian entangled state for a nonlocality test |2j and en- 
tanglement distillation Q . The photon addition is known 
to create a nonclassical state from any classical state (e.g. 
coherent and thermal states) [4| , and both of the photon- 
subtracted [B|-[8| and the photon-added squeezed states 
[H, [To| were suggested to improve the fidelity of continu- 
ous variable (CV) teleportation. The photon subtraction 
prj and the addition [12] are now practically realized in 
laboratory. 

Recently, there appeared a proposal to implement a 
coherent superposition aa) ± a) a of two product opera- 
tions, photon addition followed by subtraction (ad)) and 
photon subtraction followed by addition (a) a) [13]. Kim 
et al. particularly suggested an experimental scheme 
to prove the bosonic commutation relation, [a, a)) = 1 
jl3| , which was very recently carried out using a thermal 
state of light field [14[. A cavity-based scheme has also 
been proposed to prove the bosonic commutation rela- 
tion, exploiting interaction between three atoms and a 
cavity field [15j. Based upon the inteferometer setting 
of [13j, Fiurasek suggested an optical scheme to imple- 
ment an arbitrary polynomial of photon-number opera- 
tors, e.g. noiseless amplifier [16| and Kerr nonlinearity, 
via a combination of multiple photon subtraction and 
addition [171 ]. 

In this paper, we consider a coherent superposition of 
photonic operations at a more elementary level, that is, 
the superposition of photon subtraction and addition, 
ta + ra) , and investigate how it transforms a classical 
state to a nonclassical one. In general, a coherent su- 
perposition of two distinct operations may be created by 



erasing the "which-path" information relevant to the op- 
erations in an interferometer setting. Motivated by the 
single-photon interferometer in [13|, we propose an ex- 
perimental scheme to implement the coherent operation 
ta + ra) in an optical experiment and also show that it 
can be employed together with displacement operators to 
engineer an arbitrary quantum state in principle. As an 
example, we study the case of generating an arbitrary su- 
perposition of number states involving up to two photons. 
The superposition state Co 1 0) + C\ | 1) + Ci 12) can be used 
for quantum information processing, e.g., the nonlinear 
sign-shift (NS) gate (a basic element of the nondetermin- 
istic CNOT gate) [H, OH, and the optimal estimation 
of the loss parameter of a bosonic channel [20|. There 
have been several theoretical proposals to generate the 
superposition of |0), |1), and |2), e.g., using repeated pho- 
ton subtractions with squeezing operations [21[ , repeated 
photon additions with displacement operations [22[ , and 
quantum scissors [23|, [24|. An experimental realization 
was recently made conditioned on photodetections using 
a parametric down- converter with two auxiliary weak co- 
herent states [25]. We study our proposed scheme with 
experimental imperfections considered and show that it 
can produce the superposition state with high fidelity. 



This paper is organized as follows. In Sec. II, we study 
how the coherent operation ta + ra) transforms a classi- 
cal (coherent or thermal) state to a nonclassical one in 
phase space with the degree of nonclassicality measured 
by negative volume (area) and by nonclassical depth (26j . 
In Sec. Ill, we investigate observable nonclassical ef- 
fects, squeezing and sub-Poissonian statistics, arising due 
to the coherent operation, and propose an experimental 
scheme in Sec. IV to implement the operation ta + ra) in 
a single-photon inteferometer setting. In Sec. V, we sug- 
gest the coherent operation combined with displacement 
operation for quantum state engineering and investigate 
the generation of arbitrary superposition state involving 
up to two photons with experimental imperfections in- 
cluded. The main results of this paper are summarized 
in Sec. VI. 
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1: Wigner distribution of the state (ta + raf)\ao), where 



FIG. 2: Negative volume as a function of r on applying the 
coherent operation ta-\-rcJ (a) for a coherent-state input with 
cto = 0.01, 0.1, and 0.5 (from upper to lower curves), and (b) 
for a thermal-state input with n = 0.01, 0.1, and 0.5 (from 



+ \t\ 2 = 1, for a = 0.5 with (a) r = 1/2, (b) r = l/y/2 upper to lower curves) (|r| 2 + |t| 2 = 1) 



and (c) r = 1. (d), (e), and (f) are the contour plots cor- 
responding to (a), (b), and (c), respectively, where only the 
negative regions are colored in blue, (d), (e), and (f) exhibit 
the same size 7r/4 of the negative area in phase space. 



II. WIGNER DISTRIBUTION 

We first examine how the coherent operation ta + ra) 
transforms a classical state to a nonclassical one in view 
of the phase-space distribution. For a single-mode state 
p, its Wigner distribution W(a) is generally given by the 
Fourier transform of the characteristic function C(X) = 
Tr{p£)(A)}, where D(X) = e Xa ~ x * a is the displacement 
operator. That is, W(a) = 4, / d 2 XC(X)e aX * ~ a * x . 

(i) Coherent states: Given a coherent state |ao) as an 
initial state, the output state after the superposition op- 
eration, ~ (ta + ra^)\ao), possesses the Wigner dis- 
tribution 



W(a) = 



|to + r(2a*-aS)| 2 -|^f 



I too 



v*|2 



Wb(a), (1) 



where Wo (a) = l e - 2 \ a - a o\ 2 [ s Wigner function of 
the initial coherent state |ao). In Fig. 1 we show the 
Wigner distribution W(a = x + iy) of the output state 
(ta + ra))\ao), where r, t and ao are taken as real (|r| 2 + 
\t\ 2 = 1). In general, the negative dip of the Wigner 
distribution increases with the ratio r of the addition 
operation a + in the superposition ta + ra) . In order to 
measure the degree of nonclassicality for the output state, 
we investigate both the negative area and the negative 
volume in phase space. 



Negativity — The negative region with W(a = x + iy) < 
appears under the condition 



(x-df + y 2 < 



1 

4' 



(2) 



where C\ = \(1 — £)ao- It thus becomes a circle of 
radius 1/2 regardless of r [Fig. 1 (d), (e), and (f)], ex- 
cept the case of r = which does not change a coherent 
state (a\ao) = ao\ao)). Therefore, the negative area is 
independent of the ratio r. However, the depth of the 
negativity depends on r [Fig. 1 (a), (b), and (c)], which 
may be further quantified via the negative volume de- 
fined by V N = \(J d 2 a\W(a,a*)\ - 1) [27]. We find that 
Vn generally increases with r, as shown in Fig. 2 (a). 
For a fixed r, the negative volume Vn decreases with the 
initial amplitude ao and the extremal Vn thus appears 
at ao = (vacuum state), for which (ta + ra^) |0) = r|l) 
(one-photon state) giving Vn = ^ — 1 ~ 0.2131 regard- 
less of r. 

Alternatively, one may quantify the degree of nonclas- 
sicality by the nonclasical depth defined by C. T. Lee 
[26l | . The Glauber-P function, P(a), becomes positive 
when smoothed by a Gaussian of sufficient width, 



R(z,t) 



— [ d 2 aP(a) 

TTT J 



exp 



1 , 



(3) 



The minimum r G [0,1] required for a positive R- function 
is taken as a nonclassicality measure. Instead of calcu- 
lating the integration in Eq. (|3]) , the nonclassical depth r 
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FIG. 3: Contour plot of the Wigner distribution after the 
coherent operation td + ra) on the coherent-state input ao = 

0.5 with r = I/a/2, for (a) p^ = a) pod, (b) p |_ = dpod\ 

= a) pod) + apoa, and (d) entire distribution. 



(c) P++ + P- 



for the output state ~ (ta + ra))\ao) turns out to be 
unity regardless of r and ao by another method: if p is or- 
thogonal to a certain coherent state — i.e., (fi\p\f3) = for 
some j3 — the state is maximally nonclassical, i.e. r = 1 
[28( | . In our case, there always exists a certain /3* = — ^ao 
such that \(f3\^}\ 2 ~ \ta +rf3*\ 2 \(f3\a }\ 2 = 0. Therefore, 
the degree of nonclassicality is maximal if measured by 
the nonclassical depth r for any r ^ of the superposed 
operation ta + ra) , whereas the nonclassicality increases 
with r if measured by the negativity volume Vn- 

To gain further insight into how the coherent operation 
produces nonclassicality, we decompose the output state 
Pout into four parts, p ou t ~ V? P+- + \ t ?P-+ + r*tp. 



FIG. 4: Wigner distribution W(a = x + iy) after the coherent 
operation td + ra) on a thermal-state input with the average 
photon number n = 0.1 for (a) r = 1/2, (b) r — I/a/2, and (c) 
r = 1. (d), (e), and (f) are the contour plots corresponding to 
(a) , (b) , and (c) , respectively, where only the negative regions 
are colored in blue. 



tion after the coherent operation ta + ra) as 



W(a) 



1 + 2n 



(l + n)(|r| 2 + n) 
(l + 2n)|r| 2 +^ 



n 



+4 



(1 + 2n) 2 
n(l + n) 2 
(1 + 2n) 3 



1 + 2n 
1 



4(1 + 2 



1 + 2n 

(tra 2 +rr*a* 2 )]W th (a), (4) 



rt*p ++ , where p + _ = at Po d (photon addition), p_ + = where W th(«) = ^j+^e is the Wigner distribu- 



apoa 



(photon subtraction), p ++ = a^po^, and p = 

apoa, with po an initial state. In Fig. 3, we plot the 
Wigner function of each part for the input ao = 0.5 and 
r = -^=, where the negative region is colored in blue. It 

is well known that the photon subtraction p |_ = ap^a) 

does not create nonclassicality at all [Fig. 3 (b)], whereas 

the photon addition p^ = a) pod does [Fig. 3 (a)]. We 

see that the "off-diagonal" components of the operation, 

p ++ + p [Fig. 3 (c)], as well as the photon addition 

[Fig. 3 (a)] affects the negativity of the whole Wigner 
distribution [Fig. 3 (d)]. 

(ii) Thermal states: A thermal state p t h = (1 — 
e~^)e~^ a with the average photon number n = 
being as an initial state, we obtain the Wigner distribu- 



tion of the input thermal state. In Fig. 4, we plot the 
Wigner distribution as a function of r for n = 0.1. On one 
hand, like the case of coherent-state input, we see that 
the negative dip of the Wigner distribution increases with 
the ratio r of a) in the coherent operation. On the other 
hand, the negative region appears under the condition 



where 



C 2 x 2 + C 3 y 2 <C 4 , 



C 2 = (n + r 2 +2ntr), 
C 3 = (n + r 2 - 2ntr), 
1 



(5) 



C 4 



4(1 + n) 



2n ., . 9 

— [(1 + 2n)r 2 
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Eq. (j5j) thus describes an ellipse, as shown in Fig. 4 
(d), (e), and (f), with the size of the negative area given 
by An = Jc^ c?> • Note that the negativity emerges only 

with r above a threshold, r > yjnj (1 + 2n), and the neg- 
ative area increases with r, unlike the case of coherent- 
state input. We also find that the negative volume in- 
creases with r, as shown in Fig. 2 (b). In comparison, 
the degree of nonclassicality measured by the nonclassical 
depth turns out to be r = ^J^ 2 r using Eq. (|3j), which 
give a nonzero value for any r ^0 and also increases with 
r. 



III. OBSERVABLE NONCLASSICAL EFFECTS 

In this section, we investigate two observable nonclas- 
sical effects, quadrature squeezing and sub-Poissonian 
statistics, arising from the coherent operation ta + ro). 

First, the squeezing of a quadrature amplitude Xq = 
ae~ ld + a^e ze is characterized by (: A 2 Xq :) < 0, where :: 
denotes the normal ordering of operators. On expanding 
the terms of (: A 2 Xq :), one can minimize its value over 
the whole angle [29], which is then given by 



opt — {'■ A 2 X# :) m in 

= -2|(o t2 )-(a t ) 2 |+2(a t 6) 



2|(dt)| 2 , (6) 



and its negative value in the range of [—1,0) exhibits 
nonclassicality. Applying the superposition operation 
ta + ra) to a coherent state |ao), we obtain 



opt 



2\rf 
M 



1 



(7) 



where M = |£ao + rao| 2 + |r| 2 , and the maximal degree of 
squeezing is thus —0.25 under the condition M = 4|r| 2 . 
In addition, squeezing occurs only under the threshold 
condition M > 2\r\ 2 . 

Let a = \a \e i(f) ,t = \t\ and r = \r\e i<t>r . We obtain 
the condition |ao| = a/3|^|/ \A + 2jt|jrj cos(</> r — 2</>) for 
the optimal squeezing, which in turn gives the optimal 
|r|, given ao, as 



2 (3 + 2 v / 3|«o| + 2|c | 2 ) 



I opt 



9 + 4|a | 



(K| <>/3),(8) 



assuming cos(0 r — 20) = 1. The optimal |r| opt mono- 
tonically increases with |ao|, as shown in the plot of S^pt 
as a function of \r\ and |ao| in Fig. 5. It implies that 
the coherent operation ta + ra) achieves better squeezing 
than the mere photon addition a) on a coherent state for 
| o;o | < y/S. Moreover, it is known that the photon added 
state a^\ao) does not produce squeezing for |ao| < 1 0? 
whereas the coherent operation always yields squeezing 
[Fig. 5]. For |ao| > \/3, however, the choice of \r\ = 1, 
i.e. photon addition, becomes optimal for squeezing. On 
the other hand, a direct calculation shows that the su- 
perposition operation on a thermal state does not yield 
squeezing at all. 




|r| 



FIG. 5: Contour plot for S op t in Eq. (7) of the state (ta + 
raJ)\ao). Squeezing occurs in the region above the lowest 



Second, we consider the sub-Poissonian statistics char- 

1. 



acterized by the Mandel Q-factor, Q — — ^— j— y 
After the superposition operation on a coherent state, 
we obtain 



\ a o\ 



[M 2 +M- (l-2|r| 2 ) 2 K| 2 ] 



M[M + (M- l + 2|r| 2 )|a | 2 ] 



1, (9) 



which always takes a negative value regardless of ao,r, 
and t, with M = \ta + ra^\ 2 + |r| 2 . The negativity 
of Mandel- Q in this case increases with r and decreases 
with ao, as shown in Fig. 6 (a). On the other hand, for 
a thermal-state input, we obtain 



2n 2 (n + r 2 ) 2 -r 4 (l + n) 2 
(n + r 2 )(2n 2 + r 2 + 3nr 2 ) ' 



(10) 



from which we see that the ratio r of a) in the superpo- 
sition operation must be large enough, r 2 > \/2n 2 /[l — 
(a/2 — l)n], to observe sub-Poissonian statistics. The neg- 
ativity of Mandel-Q increases with r and decreases with 
n, as shown in Fig. 6 (b). 



IV. EXPERIMENTAL SCHEME 

In this section, we propose how the superposed opera- 
tion ta + ra) can be implemented in experiment. A key 
idea here is to erase which-path information on whether 
the implemented operation refers to the photon subtrac- 
tion a or the photon addition a\ as described below. 

First, note that the success of both the operations, a 
and a) , can be heralded by a detection of a single photon 
in an optical scheme. When an arbitrary initial state 
is injected into a beam-splitter with the other input in 
a vacuum state, the detection of a single-photon at one 
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r 

FIG. 6: Mandel Q-factor as a function of r on applying the 
coherent operation ta-\-ra) (a) for a coherent-state input with 
ao = 0.1, 0.5, and 1 (from lower to upper curves), and (b) 
for a thermal-state input with n = 0.01, 0.1, and 1/V2 (from 
lower to upper curves). 

output port heralds that a photon is subtracted from the 
initial state, due to the conservation of photon number. 
This corresponds to the action a|^), which holds well 
particularly when the transmissivity of the beam split- 
ter is large [TTj | . On the other hand, if the state \S&) is 
injected to a signal mode of a nondegenerate paramet- 
ric amplifier (NDPA) with the idler mode in a vacuum 
state, the detection of a single photon at the output idler 
mode heralds that one photon is added to the input state, 
due to the pairwise photon-creation and -annihilation 
process in the NDPA. This corresponds to the action 
a^|^), which holds well particularly when the interac- 
tion strength in the NDPA is small [12j. With these two 
schemes combined, if the which-path information on the 
detected single-photon is erased by using an additional 
beam splitter (BS2) with transmissivity t 2 [Fig- 7], the 
coherent superposition ta + ra) can be conditionally im- 
plemented. 

In Fig. 7, an arbitrary state \$) is injected into the 
parametric down- converter with small coupling strength 
K<1, which acts as 

exp(-sa t c t + sac)\i/j} a \0) c w (1 - < sa t c t )|^)a|0) c . 

(11) 

Next, the BS1 (transmissivity: t\ ~ 1) acts on the state 
as 

B ab (l-sa^c^) a \0) b \0) c 

~ (1 _ r!aSt)(l - 5 a t c t )|^)a|0) 6 |0) c , (12) 
ti 

Finally, the beam splitter BS2 (transmissivity: £2) with 
the transformations b' = t 2 b + r 2 c and c' = — r^b + t\c 



yields 

\ s \i>)) 

= [1 - ^-a(t 2 tf - rjc f ) - sa\r 2 tf + t*c f ) 
ti 

+ 5 ^aat(t 2 6t_ r *gt)( r2 gt +t *gt)]|^) a | ) 6 |0) c . 
ti 

(13) 

With the detection of single-photon at PD1 (PD2) and 
no photon at PD2 (PD1), we see from Eq. (H3J) that the 

state collapses to IV^out ~ (ta + rat)|^) a , where t ~ j^t 2 

(7^7*2) and r ~ sr 2 ( —st^). Experimental imperfections, 
e.g. nonideal photo-detection, will be further considered 
in the next section to investigate the implemented opera- 
tion particularly on generating an arbitrary superposition 
state. 

As a remark, we note the identity S^aS = a cosh s + 
ate^sinhs and S^cftS — a) cosh s + ae~^sinhs, where 
g = e §(£a t2 -<Ta 2 ) - g foe S q Ueez i n g operator (£ = se 1 ^). 
Thus, one may alternatively implement the coherent op- 
eration ta + ra) by a sequence of squeezing operation, 
photon subtraction (addition), and the inverse squeez- 
ing, &aS for the case of \t\ > \r\ (\t\ < \r\). 
However, this scheme seems to be more demanding than 
the above proposal based on a single-photon interferom- 
eter, as the number of required nonlinear resources is 
increased. Moreover, a very large squeezing s — ^ 00 
is needed to implement ta + ra) when the subtraction 
and the addition parts are comparable to each other, i.e. 
\t\ ~ |r|. Even assuming that the squeezing operations 
can be performed perfectly, we have calculated the out- 
put fidelity for superposition states Co|0) + C\ |1) + C 2 \2) 
similar to Fig. 8 and found that this alternative scheme 
does not make improvement. Therefore, in the next 
section, we will focus on the interferometric setting de- 
scribed above to implement the coherent operation. 

V. GENERATION OF A SUPERPOSITION OF 

|0), |1> AND |2> 

In this section, we show that the coherent operation 
ta + ra) together with a displacement operator D(/3) = 
e (3a^-(3 a can employed to generate an arbitrary field 
state. In particular, we investigate the generation of a 
superposition state involving up to two photons in detail. 

First, note the identity 0(P,t,r) = D^(/3)(ta + 
ra^)D(P) = (ta + ra) + ft), where ft = t/3 + rft\ which 
sequentially represents a displacement, the coherent op- 
eration, and the inverse displacement. If a single pho- 
ton state |1) is used as an input, 0(ft£,r) yields a su- 
perposition of number states as 0(ft£, r) 1 1) = t\0) + 
ft|l) + y/2r\2). In principle, a succession of 0(ft£, r), 
i.e., YliO(/3i,ti,ri)\l) will yield any desired superposi- 
tion state by properly choosing the parameters ft , U and 

n. 
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FIG. 7: Experimental scheme to implement the coherent op- 
eration ta + rcJ on an arbitrary state BSI and BS2 are 
beam splitters with transmissivitties t\ and £2, respectively. 
PD1 and PD2: photo detectors. The coherent operation is 
successfully achieved under the detection of a single-photon 
only at PD1 or PD2. 



In this paper, we focus on the generation of superposi- 
tion states Cn|0) + Ci|l) + C2I2) with practical imperfec- 
tions such as the non-unit photodetection efficiency and 
the nonideal single-photon source considered. In partic- 
ular, we investigate the effect of on-off detector that does 
not resolve photon numbers, which can be represented 
by a two-component POVM, IIo = ^ n (l — r]) n \n)(n\ (no 
click) and IIi = I— IIo (click), where 77 is the detector effi- 
ciency. On the other hand, a practically generated single- 
photon source is identified as a mixture of one-photon 
and vacuum state, p s i ng ie = + (1 - 77 s )|0)(0|, 

where r] s is the source efficiency [30]. When inserted 
into the experimental scheme of Fig. 7 following after the 
displacement operation D(/3), the state p s i ng ie is trans- 
formed to p coh = 77 S |<S| 1)/3 )(<S| 1)/3 | + (1 -Vs)\S\ ) p ) (S\o) p \, 
where \S\^) is given in Eq. ([T3]) and \n)p = D(f3)\n). 
Thus, the output state emerges after the inverse dis- 
placement D^(/3), under the condition of click at PD1 



and no-click at PD2, as p out 
Dl((3) Pcoh D a (l3) ■ ft* 



Tr b;C [p con ] 



Tr a 



z[pa 



, where p c 



IIq. Calculation shows 



where 

1$) = 

m = 
1*1) = 
1*2) = 

A = 
with R\ = 



^ (0|) 

+(1 - 77,5)77 (|*l>(*l| + ^|* 2 )(*2|) , (14) 



iiit 2 |0> + (Rit 2 P + sr 2 P*)\l) + V2sr 2 \2), 
/3*|0) + (2 + |/3| 2 )|l) + v / 2/3|2), 
(Rit 2 P + sr 2 l3*)\0) + sr 2 \l), 
(1 + |/3| 2 )|0>+/3|1>, 

5 2 |i? 1 | 2 (l- 77 + 277|£ 2 r 2 | 2 ) (15) 
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FIG. 8: Average fidelity F avg between an ideal superposi- 
tion state and an experimentally realizable state as a func- 
tion of detector efficiency 77 (on-off detector) for the single- 
photon source efficiency 77s = 0.69(diamond), 0.85(square), 
and 1 (circle), with s = Ri = 0.01. 



To evaluate the overall performance of the pro- 
posed scheme, we calculate the average fidelity of 
the produced state p out with a target state |$) = 
sin#cos</>|0) + sin sin <j>\ 1) + cos#|2), that is, F avg = 

dO sin J Q 27r d<£($|po^^"|$) over the entire range of 
angles and <p. Note that |3>) in Eq. ([T5]) is the ideal 
state that would be obtained using a perfect single- 
photon source (775 = 1) and a perfect photon- number- 
resolving detector. Thus, by identifying it with the state 
1$) = sin#cos(/>|0) + sin0sin0|l) + cos#|2), we obtain 
the state parametrization as 



*2 



p = 



V2B 2 



\Ri\ 



V2B1 



(16) 



where B\ = tan 6 cos <p and B 2 = tan sin 0. 

In Fig. 8, the average fidelity F avg is plotted as a func- 
tion of the detector efficiency 77 for various source efficien- 
cies 775 = 0.69,0.85, and 1, with s = Ri = 0.01. From 
Eq. (14), we see that even with a perfect single-photon 
source (r] s = 1) and a unit detector efficiency (77 = 1), the 
contribution of the unwanted state \<j>) to the fidelity F avg 
cannot be completely eliminated [Cf. the coefficient A in 
Eq. (|T5|)]. This is due to the use of on-off detector that 
does not distinguish photon- numbers and \<j)) arises in the 
output state p con due to the detection of higher-number 
of photons than a single-photon. In general, nevertheless, 
the fidelity F avg is largely insensitive to the detector effi- 
ciency 77, and depends substantially on the single-photon 
source efficiency r] s only. With r] s = 0.69 previously re- 
ported in [30j], a high fidelity F avg ~ 0.85 seems to be 
readily achievable within the current technology. 



VI. CONCLUSION 

In this paper, we have investigated a coherent superpo- 
sition td + rd) acting on continuous variable systems that 
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can be used for quantum state engineering. It has been 
shown that the operation creates nonclassicality out of a 
classical state with observable effects such as squeezing 
and sub-Poissonian statistics. For the coherent-state in- 
puts, the degree of nonclassicality measured by negative 
volume in phase space increases with the ratio r, while 
that measured by nonclassical depth due to C. T. Lee 
is maximal regardless of r. In particular, given an in- 
put coherent amplitude ao, the squeezing effect occurs 
optimally for r < 1, which implies that the coherent op- 
eration generally achieves better squeezing than the bare 
photon addition a) . For the thermal-state inputs, the de- 
gree of nonclassicality increases with r measured by both 
the negative volume and the nonclassical depth. 

We have also proposed an optical experimental scheme 



to implement the coherent operation ta + ra) in a single- 
photon intereference setting. Furthermore, it has been 
shown that the coherent operation combined with the 
displacement operations can be employed to generate an 
arbitrary superposition state and that a high fidelity par- 
ticularly in producing a superposition state involving up 
to two photons is achievable against experimental imper- 
fections using the currently available techniques. 



Acknowledgments 

This work is supported by the NPRP grant 1-7-7-6 
from Qatar National Research Fund. 



[1] M.S. Kim, J. Phys. B: At. Mol. Opt. Phys. 41, 133001 
(2008). 

[2] H. Nha and H. J. Carmichael, Phys. Rev. Lett. 93, 

020401 (2004); R. Garcfa-Patron, J. Fiurasek, N.J. Cerf, 

J. Wenger, R. Tualle-Brouri, and P. Grangier, Phys. Rev. 

Lett. 93, 130409 (2004). 
[3] H. Takahashi et a/., Nature Photonics 4, 178 (2010). 
[4] G. S. Agarwal and K. Tara, Phys. Rev. A 43, 492 (1991). 
[5] T. Opatrny, G. Kurizki, and D.-G. Welsch, Phys. Rev. A 

61, 032302 (2000). 
[6] P. T. Cochrane, T. C. Ralph, and G. J. Milburn, Phys. 

Rev. A 65, 062306 (2002). 
[7] S. Olivares, M. G. A. Paris, and R. Bonifacio, Phys. Rev. 

A 67, 032314 (2003). 
[8] A. Kitagawa et al, Phys. Rev. A 73, 042310 (2006). 
[9] Y. Yang and F. L. Li, Phys. Rev. A 80, 022315 (2009). 
[10] F. DellAnno et a/., Phys. Rev. A 76, 022301 (2007); F. 

Dell' Anno, S. De Siena, and F. Illuminati, Phys. Rev. A 

81, 012333 (2010). 
[11] J. Wenger, R. Tualle-Brouri, and P. Grangier, Phys. Rev. 

Lett. 92, 143601 (2004). 
[12] A. Zavatta, S. Viciani, and M. Bellini, Science 306, 660 

(2004). 

[13] M.S. Kim, H. Jeong, A. Zavatta, V. Parigi, and M. 

Bellini, Phys. Rev. Lett. 101, 260401 (2008). 
[14] A. Zavatta, V. Parigi, M.S. Kim, H. Jeong, and M. 

Bellini, Phys. Rev. Lett. 103, 140406 (2009). 
[15] J. Park, S.Y. Lee, H. Kim, and H.W. Lee, New J. Phys. 

6, 97 (2010); H.J. Kim, J. Park, and H.W. Lee, J. Opt. 

Soc. Am. B 27, 464 (2010). 
[16] G. Y. Xiang et a/., Nature Photonics 4, 316 (2010); F. 



Ferreyrol et a/., Phys. Rev. Lett. 104, 123603 (2010); 
See also P. Marek and R. Filip, Phys. Rev. A 81, 022302 
(2010). 

[17] J. Fiurasek, Phys. Rev. A 80, 053822 (2009). 
[18] E. Knill, R. Laflamme, and G.J. Milburn, Nature 409, 
46 (2001). 

[19] T.C. Ralph, A.G. White, W.J. Munro, and G.J. Milburn, 

Phys. Rev. A 65, 012314 (2001). 
[20] G. Adesso, F. DellAnno, S. De Siena, F. Illuminati, and 

L.A.M. Souza, Phys. Rev. A 79, 040305(R) (2009). 
[21] J. Fiurasek, R. Garcfa-Patron, and N.J. Cerf, Phys. Rev. 

A 72, 033822 (2005). 
[22] M. Dakna, J. Clausen, L. Knoll, and D.G. Welsch, Phys. 

Rev. A 59, 1658 (1999). 
[23] D.T. Pegg, L.S. Phillips, and S.M. Barnett, Phys. Rev. 

Lett. 81, 1604 (1998). 
[24] M. Koniorczyk, Z Kurucz, A. Gabris, and J. Janszky, 

Phys. Rev. A 62, 013802 (2000). 
[25] E. Bimbard, N. Jain, A. MacRae and A.I. Lvovsky, Na- 
ture Photonics 4, 243 (2010). 
[26] C. T. Lee, Phys. Rev. A 44, R2775 (1991). 
[27] A. Kenfack and K. Zyczkowski, J. Opt. B: Quantum 

SemiClass. Opt. 6, 396-404 (2004). 
[28] N. Liitkenhaus and S. M. Barnett, Phys. Rev. A 51, 3340 

(1995). 

[29] J. Lee, J. Kim, and H. Nha, J. Opt. Soc. Am. B 26, 1363 
(2009). 

[30] A. I. Lvovsky et al, Phys. Rev. Lett. 87, 050402 (2001); 
A. I. Lvovsky and J. Mlynek, ibid. 88, 250401 (2002). 



